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Adiabatic or slowly varying gate operations are typically required in order to remain within the qubit sub¬ 
space in an anharmonic oscillator. However significant speed ups are possible by using the two quadrature 
derivative-removal-by-adiabatic-gate(DRAG) technique where a second time derivative pulse component 
burns a spectral hole near an unwanted transition. It is shown here, that simultaneous optimization of the de¬ 
tuning and the pulse norm in addition, further reduces leakage errors and significantly improve gate fidelities. 
However, with this optimization accounting for the AC Stark shift, there is a low spectral weight pulse envelope 
regime, where DRAG is almost not needed and where the two state error fidelities are stable against pulse jitter. 
Explicit time evolution calculations are carried out in the lab frame for truncated multi-level Transmon qubit 
models obtained from a tight-binding model. 


I. INTRODUCTION 

Solid state qubits are typically formed from anharmonic 
systems such as defect centers^’^, quantum dots(QD)"^® or 
Josephson junction (JJ) devices^^^ that have unevenly spaced 
energy levels. The two eigenstates that constitute the qubit, 
are selected by resonantly coupled to an orthogonal oscillating 
field which drives Rabi oscillations between the qubit levels. 

Since these qubits are not true spin-1/2 particles, driving 
the qubits too fast or non-adiabatically leads to leakages to 
other levels in these nonlinear systems,. The adiabatic theo¬ 
rem requires that in order to remain in an instantaneous eigen¬ 
state of the Hamiltonian, the drive should be slow or adiabatic. 
This however presents a huge problem as realistic qubits are 
not closed quantum systems and are subject to decoherence. 
Hence the qubit gate operations must be performed on a time 
scale that is much faster than the energy relaxation and the 
dephasing times, while also not leaking into unwanted levels. 
The error thresholds for the gate fidelities are about lO”"* to 
10“® for physical qubits'®^'^. 

The development of techniques for the efficient trans¬ 
fer of population to unpopulated atomic or molecular levels 
is a familiar problem in NMR and laser spectroscopy'"^’^^. 
The optical Stimulated Raman Adiabatic Passage(STIRAP) 
scheme'® '^ is a well known for population transfer in a three- 
level A system. Analogous schemes has been proposed for 
adiabatic passage through a triple quantum dot system us¬ 
ing all all-electrical control'*’'®, manipulating neutral atoms 
in optical traps^" and for implementing qubit rotations^'. Typ¬ 
ically in A systems, the objective is to efficiently transfer the 
population from state |0) to |2) without populating |1). 

In the context of quantum computing, in a ladder sys¬ 
tem, efficient population transfer is required between |0) 
and |1) without populating |2) which requires adibatic gat¬ 
ing. Fast coherent qubit control has been proposed using 
resonant single-flux-quanta pulse trains^^. Somewhat re¬ 
cently, it has been suggested that high fidelity qubit rota¬ 
tions with significant gate speed ups are possible by using the 
two quadrature derivative-removal-by-adiabatic-gate(DRAG) 
technique'’^*’^"', where the second control pulse is propor¬ 
tional to the time derivative of the first. The reason why this 
works is because DRAG creates a spectral hole close to the 


unwanted transition frequency that causes leakage to the third 
level. This gating technique has been studied and adapted by 
others^^^^^. 

Consine pulses^^ and the more commonly used truncated 
Gaussian pulses'’^* have faster rise times when compared to 
an untruncated Gaussian pulses. And it is in these particular 
high spectral weight pulses where applying DRAG makes the 
most difference. There are additional known conditions to en¬ 
hance fidelities such frequency chirping'’^* to account for the 
AC Stark shift. But this is difficult to implement. Analytic 
conditions to optimize fidelity are typically based on pertu- 
bative techniques such as the Magnus expansion or average 
Hamiltonian theory. Though these are accurate to some order 
in time they might not hold true till the end of the gating pulse. 

One can ask if it gate fidelities at the threshold limit are at¬ 
tainable using simple optimization procedures and pulses that 
are just as fast, but much simpler to implement. In this paper 
it is shown here that simultaneously optimizing the detuning 
(to account for the Stark shift) and the tt pulse’s norm (to ac¬ 
count for the shift in the momentum matrix element) leads to 
much faster gate times. This has been recently shown for co¬ 
sine pulses^^. It is shown here that optimizing the detuning 
and the norm along with DRAG leads to further significant 
improvements in fidelity - but this is the case only for Gaus¬ 
sian TT pulses with larger cut offs. If one considers cutoffs that 
make the envelope more adiabatic, then the inclusion of opti¬ 
mized DRAG only makes a very small difference because of 
lower overall spectral weights. 

Another question that needs to be addressed is under what 
conditions is the commonly used three-level model sufficient, 
since the energy gaps are similar in a weakly anharmonic os¬ 
cillator. The DRAG technique can be applied to systems with 
more than one leakage transition^*. The resonance linewidth 
and the spectral weight distribution around it depends on the 
envelope. The addition of the derivative pulse envelope along 
the second quadrature, creates a spectral hole but raises spec¬ 
tral weights elsewhere. Any even(odd) parity state always 
directly couples to an odd(even) parity state in a symmetric 
potential. Moreover shorter pulses have broader linewidths 
which will more easily drive an initial state into a superposi¬ 
tion of multiple states that are nearly resonant. The increased 
chances of these parasitic transitions requires a more detailed 
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calculation of the population transfer fidelities with more lev¬ 
els. 

All calculations here, are done for the Transmon qubit. 
However the gating procedures discussed in this paper are 
relevant for any qubit formed from a nonlinear ladder spec¬ 
tra. The qubit rotations just need to be implemented on a 
time scale that is much faster than the dephasing and relax¬ 
ation times set by the decoherence mechanisms. Sources of 
decoherence for different types of qubits include hyperfine 
interactions^*^^®, phonons^'^^^, cavity mode decay^^’^^, 1// 
charge- and flux-noise noise^®^^. 

Although superconducting qubits do not have the best co¬ 
herence times, they are at the forefront for quantum com¬ 
puting because of their scalability and ease of control and 
fabrication. The Transmon is presently the SC qubit of 
choice as other SC qubits such as Flux qubits. Phase qubits 
and the Quantronium are susceptible to flux noise"*^^^ and 
the Cooper-pair-box(CPB) is susceptible to charge noise'*^’"*®. 
However CPB variants like the Transmons and Xmons are 
mostly insensitive to charge noise due to their flat bands"^^. 
The coherence times for Transmons have already approached 
the sub-ms range"^®’'^®. It is not only important to gate the 
Transmon qubits as fast as possible but ensure that they are 
stable against pulse jitter. 

In general, any practical quantum computer has to be 
fault tolerant, which requires that there be a way to either 
actively or passively detect and correct errors. Dynami¬ 
cal decoupling(DD)^°“^^ and decoherence free subspaces^'*"^® 
are examples of passive error correction. Whereas quantum 
error correcting codes(QECCs)^^ correct errors actively. 
And there are additional merits to combining both DD and 
QECCs®^ A number of recent experiments have demon¬ 
strated the working of QECCs®^ using Transmon qubits. In 
order for QECCs to work the error rate per qubit should be 
finite. It has been recently suggested that leakage errors or 
two-state-errors are the most detrimental for QECCs®* such 
as the surface code®®’^°. 

Only two-state errors are considered in this paper while dis¬ 
cussing gate fidelities. A tight-binding model is derived for 
the Transmon which gives explicit eigenvalues and eigenfunc¬ 
tions - which can be systematically truncated as required. The 
rotating wave approximation(RWA) is avoided here and the 
time-ordered unitary time evolution calculations are all done 
numerically in the lab frame. This is done for both the full 
tight binding model and the truncated models. As the detun¬ 
ings to account for the AC Stark shift get bigger as the pulses 
get shorter. One can argue that it is best to stay in the lab 
frame when dealing with large de tunings. 

The model and the method are discussed next, followed by 
a discussion of the two quadrature pulses, their power spectra 
and the main results from the unitary time evolution calcula¬ 
tions. 


II. MODEL AND METHOD 

Consider the Hamiltonian for an externally controllable 
Transmon qubit: T-L = Ho + V, where 

Ho = ^ - Ej cos{ip) (1) 

V = Eoif (2) 

where (p = 27r$/$o is the superconducting phase differ¬ 
ence across the junction, is the magnetic flux through the 
loop, = /i/2e is the magnetic flux quantum, and C is the 
junction-capacitance, Ej = Ic^ol‘2'^ is the Josephson en- 
ergy, Q is the charge operator. In the case of a current biased 
Josephson Junction, Ee = /efi^/2e, where le is the externally 
applied current. 

Since [<h, Q] = form a canonical pair and the charge op¬ 
erator Q = —ihd/d^, in the long wavelength limit H can be 
mapped on to a tight-binding model using finite differences; 

Ho = -r(4cfc+i -I- clck-i) + (2r - Uk)clck (3) 
V{t) = S{t)E, (4) 

Sx = ‘Pkclck (5) 

where Ck and c\. are the quasi-particle creation and annihila¬ 
tion operators for phase site-A:, r = Ecja^, Eq = j^C 
is the JJ’s charging energy, a = At/? is the lattice constant, 
Wfc = Ej cos{(pk) and S{t) is an externally applied time de¬ 
pendent signal. Sx is the momentum matrix. 

The calculated energy levels and the lowest energy eigen¬ 
states used to from the qubit basis are shown in fig. 1. wqi = 
El — Eq and uJi 2 = E 2 — Ei, where i?o,i ,2 are the lowest 
energy levels and A 2 = wqi — Wi 2 - For the transmon qubit 
in order to simultaneously set wgi =6 GHz and A2/woi=0.04 
we choose Ej = 22.05 GHz and Ej/Ec=l00. 

Eor easier numerical calculations, one can construct a trun¬ 
cated model by projecting out the unwanted levels as follows: 

K = (^1 Ho ( 6 ) 

r{t) = D{t)E', (7) 



FIG. 1. Potential well, calculated energy levels and respective wave 
function amplitudes for the four lowest levels for the Transmon qubit. 
Here Ej = 22.05 GHz and Ej/Ec=lQ0 results in ojoi =6 GHz and 
A2/a;oi=0.04. These calculations were done using 100 gridsites. 
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s; = (M/'|E,|vI/')=A,,(|*)(j| + |j)(*|) ( 8 ) 

where i ^ j and j > i. If only three levels are considered, 
then l^-') = {IV'o); IV'i); IV’2)} and = Xj{\j) {j + 1 | + 
\j + 1) (j|). Here \ijjj)s are column vectors corresponding to 
the state of the full tight binding model. E^ is a symmet¬ 
ric pseudo-pauli x matrix or a truncated momentum matrix. 
In a symmetric well any even(odd) parity state can directly 
couple to an odd(even) parity state. So for the three-level 
Transmon model, the ratio of the non-zero matrix elements 
are A 12 /A 01 = ^2.08. If the fourth level is included, then 
A23/A01 = 1-8 and Aqs/Aqi = 0.008. 

For the Transmon or for any similar nonlinear multilevel 
qubit, one can construct pseudo-Pauli matrices from E'^ = 
Xij |i) {j\ and E'_ = Xij \ j) (i| where j > i, therefore 

e; = -z(e;-e'_) (9) 

e; = [e;,e;]/ 2 * ( 10 ) 

The two eigenstates that constitute the qubit, can be selected 
by applying an oscillating field along 'S,x{y) which drives Rabi 
oscillations(RO) between the qubit levels and the RO fre¬ 
quency depends on the norm of the drive. E^ can be use¬ 
ful for calculating chemical shifts and modeling interacting 
multi-level systems. 

For the temporal dynamics, the wavefunctions need to be 
time evolved unitarily |^'(f')) = |^'(t)), where the 

time evolution operator is 

U{t',t)=fJ^ exp[-^H{T)]dT. (11) 

T is the time ordering operator. The time evolution calcula¬ 

tions here are carried out using a custom built C-H- program 
based on the fourth-order Runge-Kutta algorithm for explic¬ 
itly integrating the time dependent Schrodinger equation for 
unitary time evolution. Sub-picosecond time steps,df, were 
chosen so that dt l/uj where oj is the largest frequency 
component in the calculations. The resulting integration er¬ 
rors were better than or comparable to numerical precision. 
Note that all calculations are done explicitly in the lab-frame. 
If the qubit is initially in the state "ftit), then the transition 
probability to a state n at some later time is: 

( 12 ) 

Rabi oscillations obtained from time evolving the full tight- 
binding model, using 100 gridsites and a slow drive, were in 
excellent agreement with that from a truncated 3-level model. 


III. PULSE SHAPES AND POWER SPECTRA 


The state transfer probabilites can be also be understood 
in terms of the adiabatic theorem. Assume that at time t the 
wavefucntion has evolved into which is different from 
■0^(0). The leakage probability is then approximately: 




lii’njtf dH/dt 


( 13 ) 


Hence the adiabatic theorem states that a physical system re¬ 
mains in its instantaneous eigenstate (at time t) only if the time 
dependent Hamiltonian changes slowly or external perturba¬ 
tion acting on it is slow enough - which is the main design 
principle behind gating nonlinear systems. However more in¬ 
tuitive insights into the state transfer probabilities can also be 
obtained by examining the spectral weight distribution in the 
pulse’s power spectra. 

Consider a tt pulse that is resonant between levels 0 and 1 
in the lab frame 

Sx{t) = cos(cuoit)^(f) (14) 

where ^(t) is the pulse envelope. For a Gaussian pulse of unit 
norm: 


= 1 [exp(-[f-fp]V7^)-g] 

^ ^ erf{S/W)-tr.Q^ 


(15) 


Here tp is the pulse width, 7 = Wtp 12^/2 and FF is a param¬ 
eter that determines the pulse cut-off, Q = exp(— 2 /FF^) and 
erf(x) = ^ Jq exp{—t^)dt is the error function. Whereas 
for cosine pulses and shaped pulses^^’^* in general: 



(16) 


where 0 is the pulse rotation angle and ak are the fitting pa¬ 
rameters. For the usual 1 — cos, tt pulse: 9 = t :, = —1 and 

ai = 1/2. 

For regular pulses, if the pulse is too short then they be¬ 
come diabatic and this leads to leakages out of the qubit sub¬ 
space. DRAG partially fixes the problem by constructing a 
multi-dimensional pulse by adding the derivative along y in 
the RWA. The full signal will now consist of 


S{t) = Sx{t)-Syit) (17) 


where 
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FIG. 2. Gaussian pulse shapes and their derivatives in the lab frame 
for different pulse cutoffs exp(—2/IF^) 
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W = 1 W = 0.5 



FIG. 3. Fourier spectra of tp = 15 nsec pulse signals with DRAG 
(Ay = 1) and without DRAG in the lab-frame with oj = 6 GHz for 
cutoffs of {d)W — 1 and (b) W = 0.5. 



FIG. 4. Fourier spectra of DRAG pulses with a; = 6 GHz and Ay 
1 in the lab-frame for pulsewidths of (a) tp — 10 nsec and (b) tp 
30 nsec for various W cutoffs. 


here Ay is an adjustable parameter. 

In fig.2 Gaussian pulses and the derivatives are shown on 
a two axis plot for different pulse cutoffs (see Eq.l5). The 
corresponding Fourier transform of the pulses are shown in 
Fig.3 for half-DRAG (Ay = 1) and without DRAG for two 
different cutoffs of W = 1 and W = 0.5 for a fp = 15 nsec 
pulse. The drive is modulated at w = 6 GHz which is apparent 
from the plots. In the case of W = 1 the addition of half 
DRAG leads to the creation of spectral hole close to w — A 2 
(more like ^ w— I. 5 A 2 for 4% nonlinearity). Therefore when 
comparing the cases with- and without DRAG, one can expect 
an improvement in spin-flip fidelity that is proportionate to the 
net spectral dip. 

Now when considering aW = 0.5 pulse one sees that even 
though the resonance-linewidth above 10“^ is atleast twice as 
broad as that of the W = 1 pulse, the over all spectral weight 
away from resonance is smaller by four orders of magnitude. 
More over DRAG does not lead to the creation of significant 
net dip for W = 0.5. DRAG does however raise the overall 
spectral weight for both the W pulses and distributes it some¬ 
what uniformly within the frequency window. This is true 
even in the case of W = 0.5 where Sy{t) does not begin as 
abruptly(see fig.2). This obviously calls for a full quantum 
multi-level time evolution calculation in order to better under¬ 
stand its effects on the qubit subspace. 

In Fig.3 the power spectra of Aj, = 1 DRAG pulses is 
shown for pulsewidths of tp = 10 nsec and tp = 30 nsec 
for various W cutoffs. As longer pulses and smaller W pulses 
are more adiabatic, they have much smaller spectral weights 
as seen from the trends. 


IV. RESULTS AND DISCUSSION 


The gate fidelites in this section and in the next are cal¬ 
culated using the fidelity equations in the Appendix (see 
section.VII). Only two state error fidelites(F') are considered 
in this section. 


A. Three-Level System 


Consider a truncated three-level model. In fig.5, the two 
state error infidelity) 1 — F') for no-DRAG, Ay = 0, half- 
DRAG, Ay = 1 and full-DRAG = 2 is shown for the 
cases with numerical optimization and without any optimiza¬ 
tion (represented by Ay). For the numerical optimization pro¬ 
cedure the drive frequency wgi and the fundamental pulse am¬ 
plitude Ax were optimized. The frequency detunning and the 
pulse amplitude were optimized by iterating over the two pa¬ 
rameters serially. The iterations using this taxi-cab method 
were carried out till a numerical convergence of 10“8 was 
reached. These particular calculations are done for a Gaus¬ 
sian pulse with W — 1. 

It is seen in fig.5-(b), that when the detuning and the funda¬ 
mental pulse amplitude are not optimized then the half-DRAG 
results (represented by the black dashed lines) gives the best 
improvement in fidelity. This can be understood from a simple 
spectral analysis where half-DRAG burns the deepest spectral 
hole. This results in approximately a hundred fold reduction 
in the infidelity. 

However even more dramatic improvements are seen when 


A°=0 --A°=2 
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tp (units of 27 dA^ 
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FIG. 5. (a) Two state error infidelity for no-DRAG, Ay = 0, half- 
DRAG, Ay = 1 and full-DRAG Ay = 2, without any optimization 
(represented by Ay) and without optimization (represented by Ay) 
of Ax and cuoi. (b) Ratio of level-2’s population, r 2 , to the infidelity 
(1 — F'). For all calculations Gaussian pulses with W = 1 were 
considered. 
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uj and Ax are optimized (see Fig. 8 (e) and (f) for the W = 1 
parameters). It is seen here that full-DRAG gives the best 
results where the infidelity drops to 10“6 for tp = 10 nsec 
pulse. This is because the oscillating drive gives rise to AC 
Stark shift which can be accounted for by detuning w. The AC 
Stark shift energy dEs oc — 5"^ where 5 is the de¬ 

tuning and fir is the Rabi frequency. For longer pulsewidths, 
fir is smaller and hence the Stark shift is smaller and hence a 
smaller detuning is need to account for this. The shift in the 
energy levels as a function of time pretty much follows the 
signal |S'(f)p. 

Now the shift in the momentum matrix element ,(0| Yix |1) 
is inversely proportional to the shift in the energy levels and 
also follows |S'(f) p. In the absence of DRAG, the shifted ojtn 
levels that are smaller and the corresponding momentum ma¬ 
trix elements are larger due to the shifted wavefunctions. So in 
the presence of the Stark shift, Ax also needs to be increased 
in order to maintain the norm of the tt pulse (see Eq.l4). This 
is exactly what is seen in the numerically optimized parame¬ 
ters. Once the Ay is increased and the DRAG signal is mixed 
in. Ax should be decreased and the detuning, S follows the 
same trend. 

The average leakage into level-2 at the end of the pulse 
was also calculated. Qualitatively, F 2 exactly follows the in¬ 
fidelity, 1 — F'. For a quantitative assessment, it is more in¬ 
structive to examine the ratio of r 2 to 1 — F"' as shown in 
Fig.5-(b). This is shown with- and without numerical opti¬ 
mization for Gaussian pulses with W = 1. The main source 
of direct leakage is 'Pi ^2 and also virtual transitions 7 ^o-j-2 
since (0| Ex |2) = 0. The AC Stark shift induced virtual tran¬ 
sitions are suppressed by optimizing 6 and Ax- 

In order to obtain a better quantitative understanding of the 
results and to ensure that the parameter space being explored 
is relevant, it is important to compare the two-state-error fi¬ 
delities to that obtained by other authors^^. The infidelities 
for different AyS and the corresponding optimized detuning 
and fundamental pulse amplitude are shown in Fig.6 for the 
usual 1 — cos, TT pulse with 9 = n, = —1 and ai = 1/2 
in Eq. 16. These results were compared to that shown in Eig.6 
of Ref. [27] and were found to be in excellent agreement for 
the corresponding AyS (same as —2E in Ref.[27]). Overall 
the differences (in the fidelities and optimized parameters) be¬ 
tween the W = 1 Gaussian pulse and the cosine pulse are very 
small. Subsequent discussions will only focus on Gaussian 
pulses. 

The next part of this discussion focuses on arguably the 
most important parameter, W, which determines the adia- 
baticity of the pulse. The role of this cutoff parameter has 
not been discussed in much detail before. And also on as to 
what happens when the detuning and the pulse amplitude is 
optimized with- and without DRAG for a more slow rising en¬ 
velope. Two state error infidelity for different AyS are shown 
in fig.7 for W = 0.4-1.2. The corresponding optimized fre¬ 
quency detuning and pulse amplitude are shown in fig.8 for 
three different cases. 

The most striking feature here is that for W = 0.6 the 
improvements obtained with DRAG are quite minimal. This 
is because the spectral-hole burning effect of DRAG plays a 
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FIG. 6. (a) Two state error infidelity (1 — Ft) for a cosine pulse 
for different DRAG amplitudes. Ay. The corresponding optimized 
(b) detuning and (c) fundamental pulse amplitude. Ax is shown as a 
percentage change. 



FIG. 7. Two state error infidelity (1 — F') for different AyS shown 
for (a) W = 0.4, (b) W — 0.6, (c) W = 0.8, (d) W = 1 and 
(e) W = 1.2. The frequency detuning and the fundamental pulse 
amplitude. Ax, were optimized. These parameters are shown in fig.8 
for three different cases. 
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minimal role when the spectral weights are low to begin with. 
Above 20 nsec where the infidelity is flat for W = 0.6, there 
appears to be no need for DRAG. Below that limit DRAG can 
still shorten the pulse times by about 4—5 nsec. However 
in this regime the gate fidelities will be very susceptible to 
fluctuations in tp due to the sharp slope of F'. 

For W = 0.4, the behavior is qualitatively similar, although 
the pulse times now are longer for stable operations. This 
is because for W = 0.4 the resonance linewidth is broader. 
DRAG with Ay = 2.5 somewhat improves the gate times by 
~ 5 nsec, but W — 0.4 pulse is not desirable as over all gates 
times are much slower. For W > 0.6, optimized DRAG gives 
significant improvements with the most fidelity gain occurring 
for Ay = 2 — 2.5. Gate errors of about ^ 10“6 are possible 
for pulses with tp = 10 nsec. This is the fastest gate time. 
However in this regime the gates would be somewhat suscep¬ 
tible to errors in tp. 

As F' and dF'/dtp (characterizing the stability) depends 
on the region of tp, it would be more instructive to charac¬ 
terize the pulses in terms of a figure of merit(FOM) defined 
as: 

tp2 

FOM = 1 - y F\tp)dtp (19) 

ipl 

The FOM for the Gaussian pulses is shown as a function of 



FIG. 9. Figure of merit for the Gaussian pulses with various cutoff 
for tpS between (a)10 to 18 nsec and (b) 17 to 25 nsec. 


Ay for various cutoff and for pulse times between 10 to 18 
nsec and 17 to 25 nsec in fig.9. For 10 to 18 nsec, W = 1 and 
Ay = 2.5 gives the best fidelities, but this regime is sensitive 
to pulse errors. More stable operations are possible in the 17 
to 25 nsec range with W = 0.6 where DRAG might not be 
required. 


B. Multi-Level System 
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FIG. 8. Fig.7’s corresponding optimized detuning and fundamental 
pulse amplitude, A^. is shown for three different pulse cutoffs W. 


Unfortunately one of the consequences of using a small W 
Gaussian pulse is that it leads to broader resonance linewidths. 
In a weakly nonlinear system, a broader linewidth would make 
the system resonant with other transitions. Irrespective of the 
initial eigenstate, when the system is driven, the initial state 
will quickly evolve into a superposition of the states that it is 
in near resonance with. Hence it is important to go beyond 
the three-level model to see if there would be any significant 
changes in the gate fidelities. 

The calculations were redone for a four-level system and an 
eight-level system by systematically truncating the tight bind¬ 
ing model as described in sec.II. There was no difference be¬ 
tween the eight-level fidelities and the four-level fidelitites - 
the results are identical. The four-level calculations were also 
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FIG. 10. Infidelity, (1 — F') for two state errors for 8-level system 
driven by a Gaussian tt pulse for different Ay with pulse cut-offs of 
W = 0.6 and W — 1. An 8-level system and a 4-level system have 
almost identical results. 
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repeated with and without the small A 03 matrix element in 
and there was no difference as the states are far apart in en¬ 
ergy. However there are some significant differences between 
a four-level system and a three-level system which can be at¬ 
tributed to the resonance linewidth. 

The two state error infidelity is shown in fig. 10 for an eight- 
level system driven by a Gaussian tt pulse for different Ay 
with pulse cut-offs of W = 0.6 and W = 1. When com¬ 
pared to fig.7-(b) it is apparent that the addition of the fourth 
level places greater limitations on the gate times. The shoul¬ 
ders of the 1 — F' curves are shifted up by roughly ~ 5 nsec. 
This is because of additional leakages into the fourth level re¬ 
sulting from the broader resonance linewidth. For W = 1, 
the results are very similar to the three-level system results 
(see fig.7-(d)) because of the narrower linewidth. Overall, the 
W = 0.6 pulses slightly above 20 nsec would be much more 
stable against pulse jitter. This is a very reasonable gating 
time, given that the coherence times for Transmons have al¬ 
ready hit the sub-ms range^^’'*®. 


V. SUMMARY 


Simultaneously optimizing the detuning and the tt pulse’s 
norm, to account for the AC Stark shift and the momentum 
matrix elements respectively, can lead to significantly faster 
gate times by minimizing the leakage. For Gaussian pulses 
that have sharper cutoffs (larger W), optimized DRAG further 
improves the hdelity. Gate infidelties of ~ 10“^ are possible 
for 10 nsec pulses, although this is not a very stable regime. 
Much more stable gate operations are possible at tp ~ 20 nsec 
with smaller W pulses with lower overall spectral weights but 
where f™™ is more limited by the resonance linewidth. For a 
more adiabatic envelope such as this, optimized DRAG makes 
only a small difference. 

Although for smaller spectral weight pulses with fewer 
Fourier components are best for reducing leakages into un¬ 
wanted levels for an anhatmonic oscillator, for certain appli¬ 
cations such as dynamical-decoupling(DD), pulses with more 
fourier components do exceedingly when it comes to can¬ 
celling unwanted interactions and noise. For example, it was 
shown®^ that in the presence of Ising interactions and de¬ 
phasing noise for ideal spin- 1/2 qubits, second order self¬ 
refocusing pulses can improve one- and two-qubit gate fideli¬ 
ties by an additional four orders of magnitude when compared 
to Gaussian pulses. 

In order to address the wide range of issues associated with 
gating a qubit more careful pulse design considerations are 
needed. Ideally the pulses should be designed to minimize 
two state errors, minimize phase errors, minimize noise and 
be able to cancel unwanted qubit-qubit interactions. This 
would however be a challenging task and would require re¬ 
quires careful and detailed system models. 
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VII. APPENDIX: CALCULATING EIDELITIES 


More details are given here for the time-ordered time evo¬ 
lution calculations. A simple prescription is given for numer¬ 
ically calculating fidelities for all x, y and 2 : terms in the lab- 
frame for a multi-level model. As the phase terms oscillate 
with the drive, one has to go into the counter rotating frame in 
a time-ordered manner. 

Consider a truncated 3-level model. The eigenvectors are; 


| 0 ), = [ 1 , 0 , 0 ] ( 20 ) 

| 1 ). = [ 0 , 1 , 0 ] ( 21 ) 

| 2 ), = [ 0 , 0 , 1 ] ( 22 ) 

The eigenvectors for x and y can be constructed as follows; 

1±), = (10),±11)J/V2 (23) 

l±), = (|0).T*|l)J/y2 (24) 

If one starts with j'Fo) = 


{!+):„ . |-)x > \+)v > |-)y > 10)2 , |1) J, then the full fidelity at 



time (nsec) 


time (nsec) 


FIG. 11. Transitions in the lab frame between jO)^ <->■ jl)^, j—-f-t- 
and j—-f-)- |-|-)^ and the full fidelity for a Gaussian DRAG tt 
pulse applied along x. 




















































































8 






FIG. 12. Transitions in the rotating frame between |0)^ |1)^, 

|—)^ I"*"):/ l“)a; l+)x fidelity for a Gaussian 

DRAG TT pulse applied along x. 


any time t is: 


F=lTr[|Xp] (25) 

M = («'e|C/(t«)...C/(ii)C/(fo)l^o) (26) 

where the time evolution is carried out in n discretized time 
steps and U{tk) = exp{—^H{tk)dt). Here Ithe) comprises 
of column vectors of the expected wavefunctions (depending 
on the operation) and N is the number of column vectors (for 
the full fidelity J\f — 6). 

However when the fidelity is calculated this way in the lab- 
frame, it leads to rapid oscillations for x and y components as 
seen in fig. 11. An easy fix for this is to go into the counter 
rotating frame as follows: 

M = (4'e|e*^“‘"C/(f„)...e*"“‘il7(ti)e*^°‘“C/(fo)|^oX27) 

Note that that the choice between or to go into 

the counter rotating frame depends on the sign of the argument 
of U. The one with the opposite sign fixes the problem with 
the rapidly oscillating phase terms as shown in fig. 12 

For the two-state-error fidelity: |T'o) = {|0)^ , ll)^} 

Af = 2. All calculations in this paper are for the two-state- 
error fidelity denoted by F'. 


‘ F. Motzoi, J. M. Gambetta, P. Rebentrost, and F. K. Wilbelm, 
Pbys. Rev. Lett. 103, 110501 (2009). 

^ B. E. Kane, Nature 393, 133 (1998). 

^ J. Weber, W. Koebl, J. Varley, A. Janotti, B. Buckley, C. Van de 
Walle, and D. D. Awscbalom, Proceedings of tbe National 
Academy of Sciences 107, 8513 (2010). 

'* A. Imamoglu, D. D. Awscbalom, G. Burkard, D. P. DiVincenzo, 

D. Loss, M. Sberwin, and A. Small, Pbys. Rev. Lett. 83, 4204 
(1999). 

^ J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. Yacoby, 

M. D. Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gossard, 
Science 309,2180 (2005). 

® K. Hennessy, A. Badolato, M. Winger, D. Gerace, M. Atature, 
S. Guide, S. Fait, E. L. Hu, and A. Imamoglu, Nature 445, 896 
(2007). 

^ M. H. Devoret, D. Esteve, J. M. Martinis, and C. Urbina, Pbysica 
Scripta 1989, 118 (1989). 

* R. J. Scboelkopf and S. M. Girvin, Nature 451, 664 (2008). 

^ J. M. Martinis, Quantum Information Processing 8, 81 (2009), 
ISSN 1570-0755. 

E. Knill, R. Laflamme, and W. H. Zurek, Science 279, 342 (1998). 
" A. M. Steane, Pbys. Rev. A 68, 042322 (2003). 

D. G. P. Aliferis and J. Preskill, Quantum Inf. Comput. 6 (2006). 
D. Abaronov, A. Kitaev, and J. Preskill, Pbys. Rev. Lett. 96, 
050504 (2006). 

V. Letokbov and V. P. Cbebotaev, Moscow Izdatel Nauka 1 (1975). 

N. V. Vitanov, T. Halfmann, B. W. Shore, and K. Bergmann, An¬ 
nual review of physical chemistry 52, 763 (2001). 

R. Unanyan, M. Eleischhauer, B. Shore, and K. Bergmann, Optics 
Communications 155, 144 (1998). 

D. Sugny, M. Ndong, D. Lauvergnat, Y. Justum, and M. Desouter- 
Lecomte, Journal of Photochemistry and Photobiology A: Chem¬ 


istry 190, 359 (2007). 

** A. D. Greentree, J. H. Cole, A. Hamilton, and L. C. Hollenberg, 
Physical Review B 70, 235317 (2004). 

A. D. Greentree, A. R. Hamilton, and E. Green, Pbys. Rev. B 70, 
041305 (2004). 

K. Eckert, M. Lewenstein, R. Corbalan, G. Birkl, W. Ertmer, and 
J. Mompart, Physical Review A 70, 023606 (2004). 

Z. Kis and F. Renzoni, Physical Review A 65, 032318 (2002). 

R. McDermott and M. G. Vavilov, Pbys. Rev. Applied 2, 014007 
(2014). 

J. M. Gambetta, F. Motzoi, S. T. Merkel, and F. K. Wilhelm, Pbys. 
Rev. A 83, 012308 (2011). 

F. Motzoi and F. K. Wilhelm, Pbys. Rev. A 88, 062318 (2013). 

A. Poudel and M. G. Vavilov, Pbys. Rev. B 82, 144528 (2010). 

J. M. Chow, L. DiCarlo, J. M. Gambetta, F. Motzoi, L. Frunzio, 

S. M. Girvin, and R. J. Scboelkopf, Pbys. Rev. A 82, 040305 

( 2010 ). 

J. M. Martinis and M. R. Geller, Pbys. Rev. A 90, 022307 (2014). 
A. V. Khaetskii, D. Loss, and L. Glazman, Pbys. Rev. Lett. 88, 
186802 (2002). 

R. de Sousa and S. Das Sarma, Pbys. Rev. B 68, 115322 (2003). 
“ W. Yao, R.-B. Liu, and L. J. Sham, Physical Review B 74, 195301 
(2006). 

T. Yu and J. H. Eberly, Pbys. Rev. B 66, 193306 (2002). 

V. N. Golovach, A. Khaetskii, and D. Loss, Pbys. Rev. Lett. 93, 
016601 (2004). 

P. Stano and J. Fabian, Pbys. Rev. Lett. 96, 186602 (2006). 

^ T. Pellizzari, S. Gardiner, J. Cirac, and P. Zoller, Physical Review 
Letters 75, 3788 (1995). 

A. Blais, R.-S. Huang, A. Wallraff, S. M. Girvin, and R. J. 
Scboelkopf, Pbys. Rev. A 69, 062320 (2004). 

F. Wellstood, C. Urbina, and J. Clarke, Magnetics, IEEE Transac- 











9 


tions on 23, 1662 (1987), ISSN 0018-9464. 

L. M. K. Vandersypen, J. M. Elzerman, R. N. Schouten, L. H. W. 
van Beveren, R. Hanson, and L. R Kouwenhoven, Appl. Rhys. 
Lett. 85, 4394 (2004). 

K. MacLean, S. Amasha, I. R. Radu, D. M. Zumbiihl, M. A. Kast- 
ner, M. R. Hanson, and A. C. Gossard, Rhys. Rev. Lett. 98, 036802 
(2007). 

D. Taubert, M. Rioro-Ladriere, D. Schrder, D. Harbusch, A. S. 
Sachrajda, and S. Ludwig, Rhys. Rev. Lett. 100 , 176805 (2008). 

J. Clarke and F. K. Wilhelm, Nature 453, 1031 (2008). 

S. Sendelbach, D. Hover, A. Kittel, M. Miick, J. M. Martinis, and 
R. McDermott, Rhys. Rev. Lett. 100 , 227006 (2008). 

E. Raladino, M. Galperin, Y., G. Falci, and L. Altshuler, B., Rev. 
Mod. Rhys. 86, 361 (2014). 

F. Yoshihara, K. Harrabi, A. O. Niskanen, Y. Nakamura, and J. S. 
Tsai, Rhys. Rev. Lett. 97 , 167001 (2006). 

^ D. Sank, R. Barends, R. C. Bialczak, Y. Chen, J. Kelly, M. Lenan- 
der, E. Lucero, M. Mariantoni, A. Megrant, M. Neeley, et ah, 
Rhys. Rev. Lett. 109 , 067001 (2012). 

G. Ithier, E. Collin, R Joyez, R. J. Meeson, D. Vion, D. Esteve, 

F. Chiarello, A. Shnirman, Y. Makhlin, J. Schriefl, et al., Rhys. 
Rev. B 72, 134519 (2005). 

R. McDermott, Applied Superconductivity, IEEE Transactions on 
19,2 (2009), ISSN 1051-8223. 

J. Koch, T. M. Yu, J. Gambetta, A. A. Houck, D. 1. Schus¬ 
ter, J. Majer, A. Blais, M. H. Devoret, S. M. Girvin, and R. J. 
Schoelkopf, Rhys. Rev. A 76, 042319 (2007). 

M. R. Weides, J. S. Kline, M. R. Vissers, M. O. Sandberg, D. S. 
Wisbey, B. R. Johnson, T. A. Ohki, and D. R. Rappas, Applied 
Rhysics Letters 99, 262502 (2011). 

C. Rigetti, J. M. Gambetta, S. Roletto, B. L. T. Rlourde, J. M. 
Chow, A. D. Corcoles, J. A. Smolin, S. T. Merkel, J. R. Rozen, 

G. A. Keefe, et al., Rhys. Rev. B 86, 100506 (2012). 

A. G. Kofman and G. Kurizki, Rhys. Rev. Lett. 87, 270405 (2001). 


K. Khodjasteh and D. A. Lidar, Rhys. Rev. Lett. 95, 180501 
(2005). 

“ G. S. Uhrig, Rhys. Rev. Lett. 98, 100504 (2007). 

B. Lee, W. M. Witzel, and S. Das Sarma, Rhys. Rev. Lett. 100 , 
160505 (2008). 

D. Bacon, J. Kempe, D. A. Lidar, and K. B. Whaley, Rhys. Rev. 
Lett. 85 , 1758 (2000). 

E. Brion, L. H. Redersen, K. Mplmer, S. Chutia, and M. Saffman, 
Rhys. Rev. A 75 , 032328 (2007). 

D. A. Lidar, Rhys. Rev. Lett. 100, 160506 (2008). 

A. R. Calderbank and R. W. Shor, Rhys. Rev. A 54, 1098 (1996). 
5** A. M. Steane, Rhys. Rev. Lett. 77, 793 (1996). 

R. W. Shor, Rhys. Rev. A 52, R2493 (1995). 

“ E. Knill and R. Laflamme, Rhys. Rev. A 55, 900 (1997). 

A. A. Kovalev andL.R.Rryadko,Rhys.Rev.A88, 012311 (2013). 
“ A. De and L. R. Rryadko, Rhys. Rev. Lett. 110 , 070503 (2013). 

“ A. De and L. R. Rryadko, Rhys. Rev. A 89, 032332 (2014). 

^ A. De and L. R. Rryadko, arXiv preprint arXiv: 1509.01239 (2015). 
“ R. Barends, J. Kelly, A. Megrant, A. Veitia, D. Sank, E. Jeffrey, 
T. White, J. Mutus, A. Fowler, B. Campbell, et al.. Nature 508, 
500 (2014). 

J. Kelly, R. Barends, A. Fowler, A. Megrant, E. Jeffrey, T. White, 
D. Sank, J. Mutus, B. Campbell, Y. Chen, et al.. Nature 519, 66 
(2015). 

A. Corcoles, E. Magesan, S. J. Srinivasan, A. W. Cross, M. Stef¬ 
fen, J. M. Gambetta, and J. M. Chow, Nature communications 6 
(2015). 

J. Ghosh, A. G. Fowler, J. M. Martinis, and M. R. Geller, Rhysical 
Review A 88, 062329 (2013). 

® A. G. Fowler, A. C. Whiteside, and L. C. L. Hollenberg, Rhys. 
Rev. Lett. 108 , 180501 (2012). 

A. G. Fowler, M. Mariantoni, J. M. Martinis, and A. N. Cleland, 
Rhys. Rev. A 86, 032324 (2012). 

L. R. Rryadko and R. Sengupta, Rhys. Rev. A 78, 032336 (2008). 



